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COMAND - A FORTRAN PROGRAJI FOR SIMPLIFIED COMPOSITE MALYSIS AND DESIGN 

Garret N. Vanderplaats 
Ames Research Center 

INTRODUCTION 

Early evaluation of composite materials in aerospace structures requires 
an efficient means of structural sizing for a given application* It is seldom 
possible to provide simple stress limits as is customary when designing with 
conventional isotropic materials, since failure of composites is dependent 
not only on the properties and orientations of the individual plys, but on 
the nature of the loading as well. Furthermore, by taking advantage of the 
ply orthotropicity, the designer is free (within certain limits) to actually 
design the structural material through the proper choice of ply thicknesses 
and orientations, 

COMAND is one of several programs being developed in the Advanced Vehicle 
Concepts Branch of Ames Research Center to provide a general and consistent 
approach to structural analysis and design. This program is for the analysis 
and design of a multilayered composite subject to inplane loads. The principal 
method of analysis and the failure criterion considered here are those used 
by Schmit and Farshi (Ref. 1). The optimization algorithm is the method of 
feasible directions using program CONMIN, which is described in Reference 2. 
COMAND is intended to provide first level design information for membrane 
structural behavior. Another program under development includes more general 
analysis, loading conditions, and failure criterion.* 


*Program COMPOS by J. Mullen, Advanced Vehicle Concepts Branch, Ames Research 
Center . 
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The analysis and design capabilities and the basic assumptions of the 
program are presented in Section I. Section II describes the required input 
to the program and several examples of the results are presented in Section III. 
Possible future efforts in composite analysis and design are identified in 
Section IV. The principle equations used in the analysis are presented in 
Appendix A. Appendix B is a complete program listing. 

SECTION I 

ASSUMPTIONS MD RESTRICTIONS 

Program COMAND can be used to analyze a given composite panel in which 
the ply thicknesses are prescribed, or to design the ply thicknesses to 
satisfy strain and stiffness limitations. Ply orientation angles are 
prescribed, and are not design variables. Typical loading conditions and 
ply orientations are shown in Figures 1 and 2, respectively. 

The composite analysis and design is based on the following assumptions 
and restrictions. 

1, The panel is subjected to in-plane loads NX, NY and NYX only. 

Bending and out-of -plane shear loads are not considered. Multiple 
loading conditions arc considered and up to 10 independent loading 
conditions are allowed. 

2, The composite is said to fail when the longitudinal, transverse or 
shear strain in any single ply exceeds a specified limit in the 
longitudinal, transverse or shear direction, respectively. 

The composite is said to fail if the stiffness In the structural 
X, Y or XY direction is less than a specified lower limit. 


3 . 
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4. The individual ply thicknesses are designed to give minimum total 
panel thickness. Ply thicknesses are treated as continuous variables 
and several plys may be required to be of equal thickness. 

5. All plys are of the same material with the same elastic properties 
and strain limitations. Ply elastic properties (and therefore, 
those of the composite) are assumed to be the same in tension and 
compression. 

6. Ply properties are required as program input. Micromechanics analysis 
is not performed in the program, 

7. T’he composite is assumed to be symmetric about the midplane so that 
no hending-membrane coupling exists. 

8. The composite need not be balanced. That is, a ply with +45 degrees 
fiber orientation need not be balanced with another ply of -45 degrees 
orientation. Up to 18 different ply orientations are permitted, 
allowing for design of composites with ply angles at 10 degree 
intervals. Ply fiber orientation angles are prescribed and are not 
design variables. 

9. Temperature effects and temperature loading are not considered, 
except that the material properties and strain limits must be 
consistent with the design temperature. 

SECTION II 
PROGRAM INPUT 

All program input is listed heren The variables and their definitions 
are presented first, followed by data organization. No units are provided 
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for the variables. It is required that all units be consistent. That is 
(for example), if loads are in newtons and thicknesses in meters, moduli 
must be given in newtons per square meter, strains in Bieters/meter and 
stiffness in newtons /meter. 


Variables t 


TITLE (15) 
NCALC 


WPLY 

UPV 


NLC 


IPRINT 


LKK(NPLY) 


Anything may be given as a title. 

Caiculaticai control. If NCALC^O, total composite thickness 
(weight) is minimized. If NCALC.NE.O. the given composite is 
analyzed only. 

Number of plys. Up to 18 plys are allowed. 

Number of desi^ variables. This is the number of ply thicknesses 
which are allowed to change independently in the optimization 
process or the number of different thicknesses prescribed for 
analysis. l.LE.NDV.LE.NPLY 

Number of loading conditions. Up to 5 loading conditions are 
allowed . 

Print control for the optimization program, CONMIN. IPRINT = 0 
gives no print during the optimization. IPRINT « 1 to IPRII^T >= 4 
provide Increasing degrees of output during optimization, 

IPRINT “ 2 is usually desirable. 

Design variable linking. LNIC(I) gives variable number (ply 
thickness) associated with the ITH ply. For example, In a four 
ply problem (NPLY = 4), LN^ = (1, 2, 2, 3) will impose the 
requirement that plys 2 and 3 are of the same thickness. In this 


case WDV ”3. 
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X(NDV) 


VLB(NDV) 


THN(NPLY) 

EL 

ET 

GLT 

PRLT 

EPIC 

EPLT 

EPTT 

GMLT 

AllL 

A22L 

A66L 

PN(3, NLC) 


T 

Initial thickness of the design variables (IE. X = .05, .03, 

.04). If NCALC.NE.O., the composite is analyzed for ply thicknesses 
defined in X and linked according to LNK. If J = LNK(I), the 
thickness of the ITH ply is stored in X(J). 

Lower bounds on the design variables. VLB(I).GE.O, I = 1,M)V. 

It is usually desirable to set at least one VLB (I) = l.OE-10 if 
lower bounds of zero are desired, in order to prevent the optimization 
program from attempting to analyze a panel of zero thickness. If 
NCALC.NE.O, VLB (I) - 0, 1=1, NDV may be input. 

Ply orientations in degrees, referenced to the structural X-axis. 
THN{I) = Ply orientation of the ITH ply. 

Ply longitudinal modulus. 

Ply transverse modulus. 

Ply shear modulus. 

Ply major Poisson’s ratio (ply transverse Poisson's ratio, PRTL, 
is calculated internally) . 

Ply longitudinal compressive strain limit (negative number). 

Ply longitudinal tensile strain limit (positive number). 

Ply transverse tensile strain limit (positive number) . 

Ply maximum shear strain limit (positive number). 

Lower bound on composite stiffness in the structural X-DIRECTION. 

Lower bound on composite stiffness in the structural Y-DIRECTION. 

Lower bound on composite shear stiffness. 

Loads, column I corresponds to loading condition I, 1=1, NLC. 

Row J corresponds to load NX, NY and NXY for J = 1, 2 and 3, 
r^ispectively of load condition I. 
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Data Organization : 


> of Cards 

Information 

Format 

1 

Title - Anything may be given here 

15A4 

1 

NCALC, NPLY, NDV, NLC, IPRINT 

515 

1 

LNK(I) ,1=1, NPLY 

1515 

1-3 

X(I) ,1=1, NDV 

8F10.2 

1-? 

VLB (I> ,1=1, NDV (Blank card(s) if NCALC, NE.O) 

8F10.2 

1-3 

THN(I) ,I=1,NPLY 

8F10.2 

1 

EL, ET, GLT, PBLT 

4F10.2 

1 

EPLC, EPLT, EPTC, EPTT, ®ILT, AllL, A22L, A66L 

8F10.2 

NLC 

PN(J,I) ,J=1,3 (One card per loading condition) 
Begin with next set of data. Program terminates 
if 2 blank cards are read here. 

3F10.2 


This information Is duplicated in Table 1, along with a data form for 
convenient reference. 


SECTION III 
EXAMPLES 

Several examples are presented here to aid the user in making the program 
operational and to provide some insight into design using composite t«±terials. 
All examples are for a high strength graphite -epoxy composite. 

Typical ply unidirectional properties are listed in Table 2 for a fiber 
volume fraction of 0.6. The table is reproduced directly from Reference 3. 
Note that the ultimate strain limits rire not specified for longitudinal and 
transverse strain or for shear. However, reasonable values are readllj^ 
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obtained by analyzing a single ply of unit thickness, subject to a set of 
loads which are equal to the ultimate stresses. For example, given a 
longitudinal load of 180,000 lb /in. the resulting longitudinal strain will 
be ultimate strain. Therefore, a single ply composite is analyzed for the 
following load conditions: 


Load Condition 

M 

NY 

HXY 

1 

180000. 

0. 

0. 

3 

0. 

-30000. 

0. 

2 

0. 

8000. 

0, 

4 

0. 

0. 

12000. 


Note that a negative NX load is not imposed because the ultimate longitudinal 
compressive stress is the same in magnitude as the tensile stress. Therefore, 
the ultimate strains are also equal in magnitude (but opposite in sign) . 


The program input variables are now: 


TITLE: 

Determination of strains - G/E composits. 

NCALC = 1 

Analysis 

NPLY ^ 1 

One ply. 

NDV *= 1 

One thickness. 

KLC = 4 

Four load conditions. 

IPRINT ^ 0 

Not used for analysis. 

LNK(l) ^ 1 

Ply thickness - X(l) . 

X(l) =• 1.0 

Composite thickness. 

VLB(l) 0. 

Not used for analysis. 

THN(l) = 0. 

Zero degree ply orientation. 
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EL » 21,000,000 Longitudinal modulus. 

ET - 1,700,000 Transverse modulus 

GLT = 650,000 Shear modulus. 

PRLT “0.21 Major Poisson's ratio. 

EPLT « EPLC = EPTT = EPTC = GMLT = 0 - Strain limits set to ssero since they 
are not known. 

AllL » A22L - A66L - 0 Not meaningful here 
PN(I,J) - Loads, given above. 

The input data is listed in Table 3 with the corresponding output in Figure 3. 
The ultimate strains are now the actual ply strains in the direction of the 
applied load for the corresponding loading condition For example, since load 
condition 1 is the ultimate longitudinal stress, the longitudinal strain, 

EPL, under this load condition Is also ultimate. That is: 

EPLT = 0.00857 (table 2 gives 0.00870) 

Similarly, 

EPLC “ -0.00857 


EPTC = -0.0176 

EPTT « 0.00471 (table 2 gives 0.00475) 

©ILT == 0.0185 

These are now the limit strains to be used in design. 

Example 1 - Quasi-iso tropic composite 

In order to draw a comparison between graphite epoxy composites and the 
familiar aluminum materials, a simple case is first considered In which plys 
are oriented at 15 degree intervals (NPLY = 12) and subject to a single 
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unidirectional load, NX = 20,000 Ib/ln. (NY=^NXY=0) . All plys are required 
to be of the same thickness so that NDV=»1 and LNK(I)=1, I,NPLY, The 
total thickness is minimized. No minimum stiffness limits are imposed, so 
that A11L=^A22L”A66! '=0. Lower bounds on the thicknesses are arbitrarily set 
to 0.00001 in. Initial ply thickness Is prescribed as 0.05 in* The input 
data Is listed in Table 4, where the print control for the optimization 
program, CONMIN, is taken as IPRINT = 2. The program output is listed in 
Figure 4. The optimum composite thickness is 0.525 inches* The design 
is constrained by the transverse strain limit in the 90 degree direction 
(ply ntmiber 12) . The average stress in the structural X-direction (direction 
of load) in the composite is 38,000 PSl. Note that this is significantly 
less than the ultimate stress of 60, OC!? PSI for a typical aluminum alloy. 
However, the density of the composite is 0,056 lb/in.*5^3 as compared to 
0.101 lb/in,**3 for aluminum. Therefore, the relative weight of graphite 
epoxy as compared to alumintan for this example is 0.056*60000/(0.101*38000 « 
0.875 giving a 12.5 percent weight savings. 

Note that even though the 90 degree ply has failed, some additional load 
may be carried before all plys fall. Therefore, the failure stress predicted 
here may be considered analogous to the limit stress, with the ultimate stress 
being (usually) s'^mewhat higher. 

Example 2 " (0. 4-45, 90) composite design 

Due to practical considerations, ir is improbable that many different 
ply orientations will be used in most structures. In this example, the 
composite i'=s required to be balanced so that the thicknesses of the -545 and 
-45 degree plys are the same. Then there are three Independent design 
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T 

variables (NDV - 3) and the ply thickness linking vector becomes LNKT - 

T 

(1, 2, 2, 3)* The ply orientation vector is THN = (0,, 45., -45., 90.), 
A minimum stiffness of 500,000, Ib/in. is required in the structural 
X-direction. The composite is required to support the following four 
independent loading conditions: 


Load Condition 

m 

NY 

NXY 

1 

20000, 

0. 

0. 

2 

15000. 

-15000. 

5000. 

3 

-15000. 

10000. 

10000. 

4 

0. 

0. 

20000. 


The input data is listed in Table 5 and the corrersponding output 
in Fig. 5. The print control for CONMIW is set to IPRINT = 0 in this example 
and in example 3 for b7evlty. The optimum composite thickness is 0.578 inches* 
The active constraints are transverse strain limits and are Identified by 
safety factors of unity in- Fig. 5 (3 constraints are active). 

Example 3 - (0. +30, .^0, 90) composite design. 

This composite is designed subject to the same constraints and loading 
conditions as example 2, the only difference is the number of plys and 
their orientations. The composite is again required to be balanced. In 
this case, NBV = 4, IJPLY = 6, LNif *=> (1, 2, 2, 3, 3, 4), and = 0., 30., 

-30., 60., -60,, 90.). The input data and output are listed in Tablis 6 and 
Fig. 6, respectively. The optimum composite thickness is 0*532 inches and 
there are six active strain limit constraints as seen from Fig. 6. Note that 
although the number of plys and their orientations are different from example 
2, the total composite thickness is reduced by less than ten percent. 
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An additional exercise of Interest is to eliminate plys which comprise a small 
percentage of the total thickness, and solve the optimization problem again. 

For example, a composite made up of +30 and 1^0 degree ply orientations 
results in an optimum thickness of 0.526 inches. It is instructive to design 
the 12 ply composite of example 1 subject to this same set of loads, but 
allowing for different ply thicknesses (require that the composite be balanced 
for consistency with examples 2 and 3). The resulting thickness is 0.588 inches. 
Solution of this case is left as an exercise. 

Example 4 " Limit stress vs. ply thickness distribution 

In order to assess the applicability of this program to preliminary 
composite design, results obtained using COMAND are compared here with design 
curves for a (0, j45, 90) composite subjected to uniaxial tension, compression 
and shear loading (applied separately) . Figures 7-10 are reproduced from 
Reference 3. A composite with various relative ply thicknesses was analyzed 
under these separate loading conditions. Ro stiffness constraints were imposed 
a nd the lowest factor of safety was found for all strain failure criterion. 

The calculated stress was then multiplied by this factor to give the failure 
(limit) stress. The results are plotted on Figures 7-10 for 25 and 50 percent 
zero degree plys. Figure 10 compares the extensional modulus, E^. 

The results indicate reasonable comparison for compressive stress, shear 
stress and extensional modulus. However, considerable discrepancy is found 
in comparing tensile stress limits. This is because the composite is con- 
strained by transverse strain limits on the 90 degree plys. In Reference 3, 
one or more plys are allowed to fail without assuming composite failure. 
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When a single ply fails, this ply is assumed to carry no load. The composite is 
said to fail only when all plys fail individually. This again demonstrates 
the difference between the limit stress calculated here and the ultimate 
stress presented in Reference 3, The difference in results between 
these two assumptions is usually reduced when multiple sets of combined 
loadings (practical design situations) are considered, 

SECTION IV 
DISCUSSION 

A short program has been presented by viiich first estimates are readily 
obtained for design requirements of composite structures. The program is 
easily used and requires minimal execution time. Because the failure 
criterion are extremely load dependent, some judgement is n^acessary in 
choosing permissible ply orientations, so that the existence of a given 
ply orientation does not prevent attainment of an optimum design. This 
problem is much less prevalent under multiple loading conditions. However, 
it does suggest that development of an optimization algorithm capable of 
completely eliminating plys may be fruitful. 

For the results to be meaningful, it is important that this program be 
applied only to structures satisfying (at least approximately) the restrictions 
imposed in Section I, Of particular importance are the restrictions of inplane 
loading and composite synuaetry about the midplane. 

Recogni-iing the complexities of composite analysis and design as wall 
as the benefits to be gained through the use of these materials, future 
development work in this area appears warranted,** 

**Several of the topics identified here are currently being included in the 
COMPOS program by J* Mullen at ARC, 
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These efforts shotdd Include more complex loading such as bending, out of 
plane shear, and temperature loads on nonsyrametric composites. This necessarily 
requires the inclusion of more sophisiticated analysis techniques and failure 
criterion. Panel buckling under various force and displacement boundary 
conditions is also an area of interest because, with increased composite 
strengths, stiffness requirements become Increasingly important, since the 
probability of failure in this mode is increased with reduced plate 
thicknesses. Additionally, analysis and design of composites made up of 
plys of differing elastic properties is a needed and straight forward 
extension. This will provide the capability of selective re Jjif or cement of 
conventional isotropic materials as well as use of various combinations of 
advanced materials. Finally, these capabilities should be incorporated 
into a general finite element analysis and design program for application 
to large scale structures of practical Interest. 



APPENDIX A 


COMPOSITE ANALYSIS AND DESIGN EQUATIONS 


Analysis Equations 

The equations used for analysis and design are presented here. These 
equations are consistent with the assumptions listed in Section I. 

Equation numbers beginning with the letter A are consistant with Reference 1, 
The analysis is based on the ply materials properties ^LT* ^LT 

^nd , ply thicknesses, t., and orientations, 0.. 

iij i X 

The force deformation equations for the kth load condition are; 


{N}^ ■= [A] {e}j. 


[Al] 


where 



' 



1 ^xk 

11 

J-*“* 

\k 



Vk 

1 

. ^xyk 



Y 

'xyk 

L J 


0U, 
H 

3x 

9y 

3v, 

J 

9x 


3u, 

17J 


{N}, is the vector of applied in-plane loads referenced to the structural 

iC 

x-axis and is the corresponding strain state* u and v are the dis- 

placements in the coordinate x and y directions, respectively. 


NPLY 

\g ° 2^ ^^rs^^i ^ 1>2,6 [A2] 

where t^ is the thickness of the plys oriented at angle 0^ with respect to 

the structural x-axis. Coefficients (C* ), are defined in terms of 0, and 

r s i i 
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and the ply elastic constants as 


+ (C22)i nl + 4(Cgg)^ m| Hf 


[A3] 


<■=12^ = ^'=uh “i ”i “P 

+ (cjj). If m| - 4(cgg)j^ l2 m2 


[A4] 


- (‘=22^ “i "i ■" ^‘■^beh "i - “i 


[AS] 


<■=^^1 = <"iih “1 + “i 
'*' ^“^22h '*' ^^'66^1 “i *1 

(<=26h “ “i '^i + f°i2>i “i ■• ”i '^1^ 

- (C22)i ■»! 2| + 2(Cgj)^ (m. 4| - 


[A6] 


[A7] 


”i ■ ^^‘=12^ ®i 


+ ^°225i“l '^1+ Wi 


[A8] 
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where 


a. - cos 0 . - sin 0 . 

1 1 1 i 


[A9] 




"Li 


LTl'^TLi^ 




\l 1 ^Li 


\ti ^T1 


[AlO] 

[All] 




Ti 

^^‘■'^LTi’^TLi^ 


[A12] 


^^66^ i " ^LTl 


[A13] 


Note that the subscript i is not required on equations [A10]-tA13] 
since the elastic properties are assumed the same for all plys. The 
subscript is retained here for consistency. 

Given the loads the membrane strains are obtained from 

equation [Al] as 

{e}j^ = [A]-l {N}^ 

Finally the strains in the ith ply (kth load condition) are determined 
from 


lik i xk i yk i i ’^xyk 

"zik = ”i ^xk + ’^1 ®yk ~ “i’^1 Vk 

^12ik = -2=”l^i"xk + "“l^'yk + - "i> Vk 


[A14] 
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If the stresses in the ith ply are required, these may be obtained 
from the orthotropic elastic stress-strain relationships to be 

^lik “ ^^11^ i ^llk 
‘'21k “ ^lik 

''l2ik “ ^‘'66^1 ■'‘121k 

Desifin Equations 

The design objective is to minimize the total composite thickness 
(and therefore weight) ; 

NPbY 

Minimize W = S t 
i^l ^ 

Constraints on the design Include limit ply strains and lower bounds on 
stiffness. 

The limit strains imposed on the individual plys are expressed as 
constraint functions as follows: 


Gn .1 


lik ^ 

£ 0 

i = 1, 

NPIY, 

k = 1, NLC 

Ilk 


EPLC " * 






a 

^lik ^ 

£ 0 

1 = 1, 

15PLY , 

k ^ 1, KLC 

2ik 


EPLT ” ' 





G— , ■, 



£ 0 

1 “ 1, 

NPLY, 

k = 1, NLC 

31k 


EPTC 





G, 


^2ik ^ 

£ 0 

i = 1, 

NPLY, 

k = 1, NLC 

4ik 


EPTT “ ■ 





G^ 


1^12ikl 

1. £ 0 

1 = 1, 

NPLY, 

k = 1, NLC 

5ik 


GMLT 






^^12^1 ^2ik 


^*^22^1 ^21k 


[153 
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where subscript 1 denotes ply nuniber and subscript k denotes load 
condition , 

Lower bounds on stiffness are expressed as constraint functions; 

= 1. - A(a,l)/A11L <_ 0. 

G 2 = 1. « A(2,2)/A22L 5 0. 

G^ « 1. - A(3,3)/A66L £ 0, 

Constraints on strains are nonlinear functions of the design variables, 
t^. The values of these constraints are stored in vector G, (five values 
per ply, one ply after another) for each load condition In sequsnce. 
Constraints G^, G 2 and on stiffness are linear functions of the 
design variables. The values of these constraints are stored after 
constraints on strains in vector G. 

There are 5*NPLY*NLC nonlinear constraints and three linear con- 
straints on the optimization problem. Program "CONMIN" defines a 
nonlinear constraint as "active" if its value is greater than or equal 
to a specified value CT (a small negative number) , Linear constraints 
are "active" if the Jr value equals or exceeds a value of CTL. If a 
given constraint is active the analytic gradient of this constraint wa,th 
respect to the Independent design variables, t^, must be supplied. This 
information is obtained by direct differentiation of the constraint 
functions and is readily calculated using the equations of analysis. 
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APPENDIX B 

LISTING 

A complete FORTRAN listing of program **COMAND" is given here. 

In addition, program "CONMIN" is requiireu and this program is described 
in reference 2. The general program organ i’^at ion is shown in block 
diagram form in figure 11. 


COHPaSITt khALYSIS AND DESIGN PSOCRAK COHANO 


JULtj 197A 


COttPOSm ANALYSIS AND DESIGN PftOGRAH - CQHAND 


JULY, 197^ 


‘X. 



C PRLGRAh CDNANO - A PQRTRAll PRDGRAN FOR COHPQSITE ANALYSIS 10 

C AND DbSlGNt 20 

COMMON A (20l#Df (ZOIiGlSOOIrlSClSOOf f tCI EOI r At E0> 20 1 > S ( 20)* G 1 tSOO] * 30 

*02 ISOCIrC t2Q|jflSl(AOI idl20*20>* VLRC20I* WUtMZOl *SCAL(201 AO 

COPMGK YCNHNlt I Pfl ]HT * NO)/ , I TMAX *MC DN. NS 10 tr 1C «DI R*HSCAL*HFDGtf DCIt* 50 

IFOCHM* CT,CTMh*CTL*C 1 LMlH* THE TAjPHI illACrDLLFUN* DASFUN* L INDD 00 

2J>lTRMi>IU<rINF(]G 70 

CQHHON YCOnpOS/ NPLY* EL*ET*GLT>PRL)*PRTLrCPL:*EPLT*tPTC*EPTT*GniT* SO 

INLCpAllLl A22Lr Ab5L»TNtLB I* THN(IS) *CP(3)3*le]*PNC3*S l*AAI3*3l*SBI3 90 

2*ai*Ef l3l*0EPr3)*LNA(ie] lOD 

DIMENSION TlTLEtlS) 110 

EATERtiAL C0HP3 120 

C PRGGHAK FDR HULTlLATEREU COHPQSlTt PANEL UPIIMII ATIDN, 130 

C BY G. N. 9AHDERPLAATS SePT., 1973. 190 

C HASA-AMES RESEARCH CENTER, MCFFtTT FIELD* CALIF. ISO 

C RECUIPED OlPENSiaNSi TnlNPLYI * IHNINPLY >* CP13 *3* NPLV 1 * PN[3* NLC >■ IBO 

C LHKtNPLYl. OTHERS REMAIN AS NOii DIMENSIOHEO 17C 

C StORAGL kECEjIRlMENTS t DECIMAL tiORDS* CDClY ISO 

C PROGRAM - 190 

C CCPASO • 200r 200 

C CONMIN ■ o0D3 210 

C arrays * 3DC0 FOR POGRAM aS OIHtNSI'JNED IN REF. 1. 220 

C REF. 1 - VANOERPLAAIS* G. N., SCOHAND - A FGRtRAII PROGRAM FOR 230 

C SIMPLIFIED COMPOSITE ANALYSIS AND OEaImHS* NASA TM X-**,***, 290 

C AUti. 1975. 290 

t REF. 2 - SCHMIT* L. A., AND FARSHI* 9.1 JUPTIHUH LAMINATE DESIGN 20. 

C FOR STRENGTH AND STIFFNESS. INT. J. FOR NUMERICAL METHODS 270 

C IN ENGINEERING* YOL. 7* NJ. 9* PP. 519-330* 1973. 2B0 

C EQUATION NUMBERS LISTED IN THIS PROGRAM ARE FROM THE ABOVE 290 

C REFERENCE. 300 

C REF, 3 - CCNMIM - A FORTRAN PRJGRAM FOR CONSTRAlNtO FuhCTtDU 310 

C HlNlMllATlONt USERiS MANUAL* BY t. N. VANJERPLAATS, 320 

C NASA TH X-b2,2B3> AUGUST* 1973. 330 

C This PRUGFAM USES C3NMIN VERSION II* OATEU JULY* 1975. 390 

C ASSUhPtlLKS! 350 

C ailUNDART LCNOHIUNS are PRLSCRI9E0 LOADS NK, NT AND NXY. 3b0 

C ALL PLYS HAVE SAME MATERIAL PROPERTIES AND FAILURE STRAINS. 370 

C FAILURE CiclTERION ARE HAL PlY LONGITUDINAL* TRANSVERSE AND SHEAR 3G0 

C STRAINS* AND STIFFNESS LIMITS OH All. A22 AHD AQG, 390 

C VHEN ant Cine ply fails* This IS OEFIMEU AS COMPOSlTt FAILURE. AGO 

C membrane LOADS ONLY - NULTIPLc LOADING CONDITIONS ARE 910 

C CONSlOtRED. 920 

C symmetry about MIJPLAHE is ASSUMED. 930 

C CCMPCSITE need NOT DE DALAhCtO. 9A0 

C 950 

C NCALC ■ calculation LOHIRDL 950 

C 01 PC CPTtMUAllDN 9“0 

C NL.Ot D'J ANALYSIS ONLY 9U0 

C HPLY • NLPbLfi OF PITS. SYMMETRY ASSUMED. 990 

C PnOGRAM TcRMlNATED IF NPLY>0. 500 


C riLC < NUMBER Of LOAD CUNDITlGllS* 510 

C KDU > NUMBER UF ObSIGh VAKIABLBS OR INOEPENDLNT PLY THICKNESSES 520 

C FGR ANALYSIS. 530 

C IPRiNt ■ FRIlIT CONTROL FOR CDNMIN. 590 

IG READ (5*1901 TITLE 550 

READ (5*20Ci eiCAljCNNPLT*NDtf*NLC*IPRI1(T SGO 

IF (NPLY.tC.C] STOP 570 

L DESIGN VARIABLE LIRKIHG. LNMII • DcSlGN VARIABLE 5B0 

C ASSCClAItO vlTh ITH PLY, 590 

READ tbiEGCI aMKfll *1*1,NPLT1 BOO 

C INITIAL THICKNESS. BIO 

READ (5*2101 (AtlljIxItNDVI B20 

C TOTAL INITIAL LQMPDSITe THICKNcSS. B30 

D9J*G. B90 

DO 20 [■IjNPLY 650 

J'lNKdt 6B0 

UBJ>OBJ*XIJ] 670 

20 TNtn*A[j) BSD 

C LQaER bounds Cn OeSIGM VARIABLES. 690 

READ 15*2101 I VLSI II • 1*1,N9V1 700 

C PLY OfeleNTATIDM IN JECHEES. 710 

READ (V*2lCl I IHNllI rI-l*NPLVl 720 

C PLY MATEklAL PROPERIIES. 730 

HEAD (6*210) lL*ET*GLT.PRLT 7Q0 TO 

PRTLoPPtlAETyEL 750 O 

C STRAIN AML STtvFMbSS LIMITS. 760 

READ I5,2l0l CPLC*bPLT*EPTC*EPTT,GMLT*AlU*A22L*A66L 770 

C LOADS FOR EACH IJAD CONDITION. 7BU 

DO 30 I«t*NLC 790 

C LOADS - NX* NY AND HXY FOR THIS L3AC CDNOITION. BOO 

10 RLAU [L*2IUI [PNU»n*J-l*3I BID 

C NCON > KUMBLR OF CDNSTRAlHTS. B20 

NCCN»5*NPLY*NIC 030 

DC QO l*l*NCC(i B90 

All isc(n*o BSD 

NI-HCCM*! B6Q 

lSC(hll*l 070 

lSCtNl*lJ-i 060 

lSCtNl*2)*l 090 

If I AllL.r-T.l.OE-lO) NCON*, (CON* 1 900 

IF [*''cL.GT. 1.0C-10I NCJN*HCiINtl 9X0 

IF (.vctL.GT.l.OE-lOl NCUN*NCDNtl 920 

C PIIHT input INFUfiHATltH, 930 

IF IIiClLC .EC.t» RRITEl6(A6flr 990 

IFCHCAIC.M ,( 1 RRITEIF,a70I 950 

vRilE (G,22CJ 960 

vRITt (6*1501 TITLE 970 

HBITL (6*2=0 HPLT,NLC 960 

VPJTE l6i2AO EL*ET.GLT*PRIT*PMTL 990 

VRITE (6*250 loco 


CETHPG^lTt. ANALYSIS Attti UESIGN PMJCRAti - CilHANO 


JULY I 197^ 


LUHFaSlTb ANALYSIS AHU OESIbh PRQGRA'1 ~ CQHAND 


JULYf 1<»7« 


Oil !>0 I>ifhPLY 1010 

5C HRIYE ffr»770l I*TN(M,THrnn»lNK[n 1070 

GHLU*-iihLY 103C 

bAITb larZoO) tPLCrEPlTiEPTCrirprT.OnLTljGnLT 1040 

liRnE (L»Z90) AlXLf A72L« AtiEl lOfril 

HRITE lb»300t lOEO 

00 bO I-ifME 1070 

tiO MRlTt (Ef3I0> I«(PNU>I>>J>1«J) lOOu 

C 1M11IAUZE CONtU;! PARANbTEAS TJ OEfAULt VALJE^. 1090 

llnAKKlO 1100 

KSlOc-1 1110 

IChDl»«0 1120 

bSCAL'O 1130 

nFUG-5 11^0 

FDCh«C. 1150 

fOLHn^Ct 11*0 

CT.C. 1170 

CTniN>C. 1150 

CTL«Q* ll9fl 

CTU*Ih>0. 130C 

THtTA-C. 1210 

PHI>G. 1270 

DeLF0N<0* 12 SO 

DAEfON'O. 1240 

IZSO 

lTSti»C 12tO 

C CONlftiiT PLY angles to RADtANi. 1270 

Uo 7C l'l»r.PlY 12B0 

TKHII 1-TrMlS/:(7.2»ti?76 1290 

C UPPER bOUNLS ON OcSlGK VARtAOLcS AROITRAPILY aET • LOO • 1300 

70 WUBUI«100> 1310 

C PLY STIFbr.biS CuuflClENTS. 1320 

CALL ClMPl (NFLY»THN*flf£ r.GLI»PRtT»PRTLf CPI 1330 

AN1*2C 1340 

hN2*5« 1350 

Wi3»20 1360 

NHA«7t 1370 

SH^>4C I36u 

IF <NCALC.bC. 01 CALL CjnNIn IC JnP3i <13 J r OF t Gj I SC r 1 C • A j S j G1 >G2i C*M 1390 

*SljdrRLEf VLBfbuAL iNNl iNS2iNN3<Nll4f NN7I lAOD 

C PKIM ANALYSIS REaULIS. lAlO 

bRITE IbfSitrl 1920 

I.R1TE IbAiSGI tlfLt 1430 

VRITb tc>f330l 1440 

P?*£CC.^ubJ 14S0 

C PLY THlCAKLSStJ A9U PEKCt'iT uf TOTAL ThUKNEsS. 1460 

Oa aO 1»1iNHLT 1470 

J-Lhn<ll 14H0 

TNtn»*lJI 1470 

PCt«PP*Thm 1500 


AO kRITt IbjZbll IflN<n>PC1 1510 

C total ItiiCKALSi. 1520 

kRlTb tbj34Gl ODJ 1530 

C SriFFhLSBLA AND FLt LASl L I T lES . 1540 

CALI C0NP2 INPlTj T ritCPt AA.B91 1590 

C PLY SIRaIn* Ah£. C3,1P3SlIb STRESSES FDR ALL LOAD CONDLTHmS. 1560 

fcRJTE <6r35U 1570 

OD JOL Ll>lfNLC 1500 

kRlTE (6»36CI II 1500 

C COPPBSITL STRAINS- 1600 

EP (lT°<ld(lfl)*PNt jLfm«>9t!(lr2)*'PN<Zf II >>BB(1 a3I*PN(3(I1I 1610 

EP<2l>bBl2fl]*PNU»Il)AdB12»2}«PrU2f IllAEi[l(2r3>4PH(3(UI 1620 

LP 13l«QDt3tl)4PStl,II M'BStlAZIAFtllZf m*0B(3»3l*PH(3(ll) 1630 

C PLY SthAINS ANU SEFcTY FACtURS. 1640 

CJ 9C J<ri<hFlt 1650 

SnbTA»TH(ltJI 1660 

AL«C0SITriLTA1 1670 

Ah«5lMTnblA> 1600 

ALZ»AL*AL 1600 

ah2«AK*ah 1700 

C STSAIKS. 1710 

EPl«ALZ*Lp<Jl»Ar1i*ePIZI*AL*AhAtf < JJ 1720 

EP2*AhZ*bP<ll*ALZ*tP<2l-Al*AJ<*LPt3t 1730 

bP3»2.*ALAAf5*<tP(Z)-EP< 1H*IAL2-A!12 JAfcPIll 174<J 

C SET SlFAlNS TO NlhlflliK ABSOLUTE VALUE OF i.JE-20 TO PKEVEHT 1750 

C CIVlOt PT IbUL'. 1760 

IF IABS<IP1>.LT-;,Ue-2k» EPl«1.06-2t 1770 

IF lAUSltPZS .L 1 .1 .0 e-2l,I £P2«1.0fc-2O 1700 

IF tAUStLP3I.Lt.l.GE-20] i?3«l,0E-2C 1790 

C SAFLTT FACILM- 1000 

SH>£PLC/lPi lOlO 

SFZ«LPTC7tP2 1820 

SF3»-GMLlibP3 1830 

IF tEPl.Gl.0.1 3Fl»EFlt/EPl 1890 

IF (EP2.CT.0.) SfZ«ePTT/EP2 1050 

IF <bPj.CT.Gi> SF3«“SF3 I860 

IF ISFl.uT.lCC.I SFI-1C3. 1870 

IF <SF2.i>T.100.l SFZ-100. 18B0 

IF tSFi.GT.lOt,) iH-lQO. IB90 

WPITt I6.37CI J*ePltSFi,bPZiiSf2jEP3»SH 1900 

90 C3NIINLL IHO 

c cuftptism sifcAiNS. 19Z0 

kRITt <Of3tCl < EP<K| ,A«1*31 1930 

IJG CQMINLfc 1940 

C COHPQSHE STPtSSEj, I960 

kPITb 10.39&S I960 

CC 12L 1*1, ML 1970 

UU lit J"li3 1980 

110 CUl-FMJill/U'jJ W90 

kRITE 16,4001 It <G<J ), 2000 


CJCPCSUE iMitrSli 4KD U6SIGN PRf.RAH - CDHmO JtlLY» 1S7S 

120 CONIINUb 2010 

C COMPOS] rt STItmtSSii* J{j 2 y 

WRITE I6j160] 2030 

$f*|.00, 20W0 

IF (AllL.GT.lr 5P>A«(l»U/Aj,U 2030 

If (SF.or.loCiJ Sf>10L» ZOtO 

VSITt lb>I7Cl AAUiltiAUL*SF 2070 

SF-100, 2000 

XF (A22L.G1,l»0t'l0) SF>AA(2f2I/A22L 2000 

XF (Sf.Gl.loO.) jP*10C. 2100 

h»m (btlbC) AAI2>?]rA22LpSF 2110 

3F*10C> 2120 

If tAttL.GT.X.OE-UiJ SF«AAt3f 31 XA6bl 2130 

Xf (Sf.Ul.lCO.l SF*tOO. 21W0 

WftllE (b>1901 AA43pi]>Aa6L«SF 21^0 

nRXTt (biAlCil 2lbU 

DO 130 I'ijj 2170 

DO 13C J<1j3 2100 

as U» Jl»l]BJ*EU ( i*J) 2100 

130 AAUr J I'AAIXp J WJ3J 2200 

C CGRPOSlTt StkbSl-STRAlN PbLATijMSHIPS- 2210 

WRXIl (&jA2u 1 HAAtii J1 >J«Ip31>I>1p3) 2220 

WRXTfc (bj'>3Gl 2230 

C COhPUSITE ilftilN-iTAfiSi RtLATxONSHIFS. 22^0 

WRlTk (bjWAGI ((3UIIp3ljA«lp3Wl«lt31 22D0 

t COHPDSIU ELASIIC COMSIAHIS. 22t0 

IXiI./batlill 2270 

tT^i.yftnr ?#i) 22ao 

(iXT*lt/BlH3j31 2200 

PrlXY*-bd l4./2|ybB( L> L ) 2300 

PRTA«-Ub u,2) Xu(U2p21 2310 

write LApEYiOXTiPRXYpPitTX 2320 

CJ TC 10 2330 

C 23A0 

HO POkhAT UiAxl 2350 

IbO fUkhAt C yiWAp ^ mTI I tbX 1WX( IjaAl 23b0 

IbG fUHPAl C//2>Xji0rtL.0)1fDiIT£ ftlPiaAiML ST 1 f FNl S 3cS/ 27X r ottAC TUAL >7 Xp OH 2370 

lREU.IDLU/i7XjSHVllUi;>7X>9HVAL0c>t)X)^HSiFit 2300 

170 F0RJ1AT ll4Xp3h4UttU.3»lAabl3.Sp3XpF7.2l 2300 

lot FOSrAT U*f»p3nAZ2pil3*»,lA,E13»i» J*pF7.21 2400 

HO FOkKAT UWXp3>1AbapEX3.:rplXpL]3pbp3XpF7.?l 2410 

200 FOBPaT Iltlil 2420 

210 FOPPAT tiFiC.21 2430 

220 FOAAAH 7/41Xp2HGF/y )CXp23HSTI!n2 TAU COHPOSXrE PAPELl 2440 

230 FflAfiAT OF PLTj I5p /p 27Xp 2!iNND. OF LOAD 2450 

1 CChblTIL^5 >pX5I 2400 

240 fOSpAl t //25 Xp35HPLy PKbPEATItS - ALL PlYS lOt MT ICAL / 26X» 22HLDHCIT 2470 

lUOlbAL PULULU5 •pLi2.5/26Xp22HTRAN3v£A5tl MtlOULlli «p £ 12 . 57 26X p 2 2H 24U0 

2SHtAp FODliLllS •pU2.5/2bXp22HPaiJSOr4<3 AAIIJ, l-T •»EX2.5/2 2400 

3bXp22bPai55LMS RATJUp T-L •p£12.51 2500 


V 


250 

260 

270 

ZliD 

2'JO 

300 

310 

320 

330 

34 c 

350 

36C 

37L 

3bO 

344 

400 

410 

420 

430 

440 

450 

4bc 

470 


CanPOSlTL AhALTlIS AtiO DESIGN PKOGRAit > CDhANO 


JULTa 1974 


FOPCAl (y p 14Xp3lriPlY THlCKNtSStSp 0HXbv*TAtlQNSyp27HANU DESIGN VARl 
lABlt Min3£KS/20Xp47HPLT NO. THIOHbcSS THETA DES« VAR,. NO 

FJPNAI (21Ap:4p3XpE13.5pF 10p2) 

FORRAt U5Xpl9p3Aptl3.5p3XpF7.2p5Xpl5l 

FOKRaT iyy3H,l7riPLY STKAIN LIrtIIjyl4Xp23.UONGiruOINAL STRAIN.GEpp 
U13.Sp4H .ANI,,Lc.pE13,5714Xj 23H TRANSVERSE STRA IN.uE, pbl3« 5>9H .A 
2hO*Lcp ptl3.5yl4Xp23H SHEAR STAAIN.GE .p£ 13. 5p9H . ANt..LEp>E13* 

351 

FJAHAl Wy^SXpibHSnFFNtSS L IN ITS/ 33X p THAI I .Gb ,p £12 . S/ 33X*7HA22.CE 
1.#£12.5/33 Xp7hA65.0cppc12.5I 

FUkHAT [/y40Xp5)IL0ADS /l5Xp IOHLOAO CONO. pLXpZHHXp 12X» >''*iYpl2Xt 3KHVV 
i> 

FiiPRAl il5x#ltp6Xpfcl2.5p2Xpfci2.5*2Xp£12,5J 
FGAtiAT UHIpZoXp jOHOt-SX-pN ANU/Ok ANALTSIS RESlilTSl 
fOARAT ty/yiSXalSNPLY INHQAHATI0n/25Xp32HPLT NO. THICKNESS 
IPtRCbMJ 

FUPhAT I35 Xp1lH — plXpluM — /23Xp HhTlIIC KHbSS -pElZ. 

1 5p 4XpCHI CC .CO ] 

fORRAl |yyyy37XpllPtPLY STKAXNS/33X,23rt5,F, • SAFETY FACT0R/33 Xp 26B 
lEPL * LIlfc&ITUOlNAl STAAXH/33Xp24N£PT * TRANSVERSt StBAIN/33X»19H 
2EPt( > SHEAR SIHAIH) 

FtlRP-AT ( //ScXpIOHLOaO CullD.p X5/5Xp7HPlY 40ppbXp3h£PLpbXt4NSpF«peXj 
l3HbFT/bXp4HS,F .pdXpAPitPLTp 7Xp4HS.F. > 

FOWhAT t 5XplS,lXptl2.5p2xpF7.3.2XpE12,5p2XpF7»3p2X,bl2.5p2X,F7,3| 
FOdRAT [ y tXp47riL0PiPaSl IE SIKAHS AbFERbKCLO TO STRUCTURAL AX|S/SXp 
ItrilFX • pE 12f : plXpbHt FY ■ il l2.5p3Xp7-<EPXr » rbl2.5] 

FOWNaT I //2iXp30HHtHBRAr,E SIRcSSES IN CunPOSX tE/i5Xp IOHLOAO COHOp* 
15Xp7hSlGnA-XpbX. 7H3 XGHA-Yp 7XpbHTAo-xYJ 
F0MPA1 tl7Xpl5pTXi3£:3.5J 

FuhFAI I /y2tXp42riCOtP-ICXENTS ill- STAF iS-lTRAlN R£LATlaKSHlP5/2BX,26 
IMBILAUB IG STkliCrUHAL Axb>) 

FDWBAl llOXptHCll - >L12.5p4KjoHC 12 ■ P tl2.5p4X, 6ltCl6 • p£U.Sp4Xp 
l/32>p£FC22 • p : 12 .Sp 4Xp5HC2b ' Ptl2. 5/20KpVHSYHHbTitICp25Xi BHCbb ■ 
2pt;12.5} 

FOkHAT 1 yy£UAp42rtLDbt 1C.6NIS JF STAAlH-SiKbSS RE L AUDNSKIP 5/2 BXp26 
IHALLATLO TC STPUC10AAL AXES) 

FJRHAT liOpEHail . pbl2.5p4XpcHQ12 • p £ 12 . 5p 4Xp bhjlb • «E12.5f4Xp 
l/32XpChQ?2 • . tl2pSp4Xp(,h02b - p b 1 2 . 3/20 X pVHSYHNETRXCp 2SX«6HQ£6 • 
2pE 12.51 

POPHAI (y/ZbXppTriCOKPDSlTfc cLAiTIC CU.vSTAHli/llX . VMfcX « pElZ.SpSXp 
l5HtT • n.l2,5i4XpbH3XY » ptl2.V/9Xp 7HNUXT • ■ E 12. 5p 3X p 7HNUY X » pEl 
•2.5) 

FO AHA I ( it,; p 3<tx p bHOtS lOH 1 
f 0A/lAmHl,37Xp3HANALTSIS I 
End 


25X0 
2520 
2530 
2540 
2550 
2560 
2570 
2500 
2590 
2600 
2610 
2620 
263C 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
2770 
27B0 
2790 
2BOO 
2b lu 
2020 
2830 
2640 
2550 
2660 
2870 
2BBD 
2590 
2900 
2910 
2920 
5930 
2940 
2950 
2960 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 



SUHPOLiUNt tbhPi (hPlV> THN>lLrbT(SLr>PHLrrPRrt.jiCP] 

01PE«il(iN CP(j>3^1«»> iHMiaJ 

RDVIINE IC CAUCt^LArt PIT SIIfFNESi COEFICIENIS - Alt PLrS THE SAME 
tLASlIC PRUPtfcItES. 

OY b. N. 9AM)EPPLAATS sept.* 1973« 

hASA-'AnbS hbSEARCtl CcMTi-P. HOFFcir FlEUJt CaUF. 
frATlRlAL elastic i-ejPlRI'ES. 

Pftl«l,/(1»''?PLT*PBTL I 
taUATICa Alo 
C11-U*PH1 
Ei}IjAT}CN All 
Cl£»PklL*tl.*PRl 
LJLATtQN AlZ 
C2t'ti*esi 
FOLATIUI A13 
Cbt>«GlT 

CO FUR ALL PITS, 

CO 1C l»l*KPLY 
INl TA>THM II 
tmiATlliN A9 
AL>CaS(TriElA3 
AhiSIMIrlETA) 

AL2«AL*AL 

At3-AL»AL2 

AL4-ALi*ALJ 

An2»ApAAM 

A43>A«**At 

AH4"AP?»A«i 

tnUATlLTa A3 

CPtl»l*n'Cll*AL4 + 2. *ClZ»ALi;*4f1Z*;22*An4A4.*C6o*AL2*Al1E 
taUATUH AA 

CP ^•ClJ•Al^»A^^♦C^^•^ AL4 *AK<i ^♦C^^•Al^•4(1^-••,«Co6*AL^♦AM^ 

EOUATICn as 

CPIlrirl J - tll»AL3«AfHCt2*lAL*Art3-AL3*AS»-CA2AAL*Allj+Z**C66*<At»AP3 

i-an«AL3i 
‘OUATILN At 

CPiZiJ* J )*CU*An4AJ.*ti£*AL2»An2+:Z2*AL4»4.*Ct6A*LZ*A«Z 
eOLAllLN A7 

CP I2t3f [ l*Cll»AL*4A3tCl Z* tAL3*Att-Al»AM3l-CZ2*AlJ*Ant2,AC66A I At 3*AM 
L-AL*AF3I 
taLpATlCA Ae 

CP I3>3r ^|•CU*At^•A,■^^-^.♦Cl^»ALZ*At^^>t^Z*AtZ•A^^♦CbO*^ Al4-2i*At2*4 
iaZ»AH41 

IHPOSt iYpMIBY ON C». 

CPlZilrl l>CP{l>ZTn 

CPI3f l«l tk^Pllf3>ll 

CPI3>Z»I)>CP1Z>3|1 I 

CONI INU 

KtlURS 

£N0 


10 

20 

30 

40 

50 

60 

70 

flO 

90 

100 

110 

120 

130 

140 

150 

160 

170 

ISO 

190 

200 

210 

220 

Z3t> 

240 

250 

260 

Z70 

260 

29g 

300 

310 

320 

330 

340 

350 

360 

370 

3B0 

390 

400 

410 

420 

430 

440 

450 

460 

470 

4B0 

490 

>00 
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SOBRDOIINF C0hP2 [NPtY»TN*Cf f A>ei 10 

OINENSIJN A[3r3lf 0I3 p 31> TMl2l* CPI3 a3»1BI 20 

C ROOTINI TD CAi.CULATe KEHaftANt SlIFFNESStS ANO FttAlBlLlTlES OF 30 

C COrtpOsm NAtiE OP OF NPLY PLTSj EACH «ITM THE SAKE HATEHIAL 40 

C PROPERTIES. SO 

C BY 0. h. VANOERPLAATS SEPT.j 1973, 60 

C NASa-apES RtSEARCH CENTER* flijFFETT FIELD* CALIF. fO 

C STIFFNkSS COEFICIENTS. BO 

C ZERC A. OO 

DO 10 1>1*3 100 

LJ 1C 4*1*3 110 

IL AI1*JI*0. 120 

C bJlLD A BY SL'Pb*<POSlTlQN. 130 

C LQLaTIGN A2 140 

CO 30 1»1,NPLY ISO 

T.TNin 160 

00 2C 170 

CO 20 H.J.3 IBO 

20 Al J*Nl<Al j,iil*T*CPtJ*KrII l^D 

30 COhlINLt 200 

C IHPDSI SYnitbTMT 34 A. 210 

A[Z.iJ«All«21 220 

A{3*1I*A(1>3I 230 

Af 3*21*412.31 240 

C FlLAABIlITY CDeFICIENTS - INUtHit OF STIFFNESS. 260 

C BUILD E>a-1nvLnS5. 26U 

Dc 1>A[l*l}*Al2.2l«AI3.3l*2.*A{l*2l*At l*3]*A[Z,3l''A|l.l]<A(2i 31 *AI2 270 

E, 3 l-AtirZJ *AI 1,3 l*Al 1,3J- A( 3*1I*4U,21*AI 1,21 200 

□ET»l./LlT 2^0 

B( I* 1 1*021 *(AI 2* 21 *AI 3* 31-4{2,3 1*412.31 1 300 

Bl l,2J-0LT*|AU,3|AA<2.3>-Ari,ZI*Al3,3l I 31u 

Bl 1,31 >ut 1*1 All* 2)«AI2.3)-A< 1,31*41 2,2 M 320 

b:2,2I*OU*IA1 1, 1]*AI3,3]'>A<1,3|A41 1,31 1 330 

9l2,3><OLT*lAlt,2)*Atl,3l-Ali,lJ*AI2,31 I 340 

01 3,31 *BcT*IAll,il*AIZ,2l-4(i,21*A(l*2I I 350 

C IHPCiSt SYPHtTht ilN 3. 360 

£I2,1)>BI1,2I 370 

at3,l)>Ba.3l 360 

B(3,2>>B|2,31 390 

RcILiiN 400 

t»D 410 


ro 

1.0 


c 

c 

c 

c 


COhPCSirt *hALtSlS iuD DLStCH PMOGIUn - lCKJIMI} - COnP3 JULY* 1S7^ 

3UeSCbTtHf CLKPJ <lHf t,*UBJ.**OF»G.IC»Af 10 
Ca/'HON /CNPHl^ IPC trtfjNDVf ITH 4 **«C 0 K*MS laL# lC^O^R•MSCi^L^NFa&,FOCH^ 20 

LFOChHfCFf CIRliirCTLfCTLhXNi THL YAl r PHI * N3C * P£LF U^l* DAUFUHt L INO 30 

2BJ,iTttU,lTE:P>lHFaG 1,0 

DltlLNSION KtM »*DF (HI l,GtM2l > 1C (■!31 1 A |Ji3i M ] 5 q 

C^ItRhAL P0UT1HL FUJ CONBIH FOR CJI«Hait7!, FANtL CLSICtt. 60 

BY G< U, BAhCtPPLAAII SEPIx 1973. 70 

WASA-APtS fctiLARCH CtHTbB, HCFFtll UELOi CALIF. SO 

This II A bLFFcR BETitLEN CDNMN AriO COnPA. 90 

CALL CL6PA ( tPiFa<OBJ*ltDV>CT*CTL>r<AC« KiOFi G. A* JCir<I*N2*N3) 100 


/ 


/ 



C 

c 

c 

t 

c 

c 


10 


2C 


JC 


AO 

^0 

c 


c 


c 


CuBPOS. U AhALYSIS ABD OESI&S PHOSBaM - COHANQ - CCHPA JULY* 197A 


1-iOFULItr.E CChpa I 1.HF(3*QD JiHdV.C F.CILjNAC* i(iOF*<;*A»iC*Nra*N 
♦L. .HH3 I 

CORHiIh FCGRPtS/ nPLY»FL»ET*GLT,PRLT»PRTLniPLC*EPLT.ePTC«EPTtiGnLT* 
lALCrAlU»AC2LiA66LtThrifal*lKhll6liCPt3r3.ie)iPHU*ai>AA(3>3]*Bai3 
2*3I*LPI3Ji[ilPM]iLrKlltf I 

DlnFNllJa InPOJ. KIH)IU,aF(NMI*G(NN^IfK(NN3]f A(HII3*NHlt 
ROUTINl 1(1 CALCULATE FuriCiau VALUE* COnSTRAlNl VALUES AND 
(.AAUUAr U lUNCnDH AAU ACTIVE CONirriAIltTS F£lS CCnPOS lU 
AIIALYSU Arc DESIG .1 PFOGRAK - COnAMP. 

BY G. A. VANOERPLAAII ScPT.t 1973. 

K4SA-AMi FtitARCH CthTEfi, VCHHETT FIElB* CALIF. 

IF (1GFU.GI.2I CO rO 3U 

L3JLCTIVE 

uB J *C . 

(JO 10 I*1*NFIY 
J*L6KII) 

T>i(l)*t(j] 

CBJ«LfaJ*I6(l) 

IF llMU.tO.3l RtTUV*! 

IF lihFu.U.2] CJ rj BO 
C3MIAU 

GIaOIlM CF GBJFCTlVi 
£0 3C I*1 *NCb 
t'tlJpt, 

CO BO t«l*BPLY 
J-LSKt n 
OFtJl«LF(j|ti. 

IF ( lAf Lj.tL.3J RETLRS 
COFlUNLL 

COMTRAtKTS ARu OHADIEhr OF ACHVE C TNil-t A (NTS . 
f.CTLl*C 

Ir IIt.F J.LCi B] HAC«0 

SriFFNESS AM FLcAXd [L (T (ES. 

Call Cl'BPP IhPLTrTHi LPi AA.3U I 

CD 17L J-iiALC 

IKViRIt OF iLUAnUt Al. 

£P< 1 1*bSU*ll *P«U*J iFoMl l,Z 1»PN( ?* I »*SBI U J |*Put3*n 

IP (2l-btp(L) 11»PM J *n*BB< J*2 )*PN( J* 11 * 00(21 31 *P>l(3rI] 

Ek(JF*bBIj* 1]»PHI t,l HillH 3*2F*P*J(2, n«BQl3, 31 *P 1 I (3,n 

□0 17U J'lfhPLY 

FHClA*1hi<1Jl 

AL*CCS(TBt 1 Al 

AM.Slf-dBETA) 

AL2'A1*AL 

AH2«AB*AP 

E.UAIICM AIL 

EPl-AtF»lPm+Ar?»i:PlJJ*AU*AP*Ep(3F 
IP2»AP7*E PU J*Al2*tP( 21 “AL*Ah»cP[ 3 1 
FP3-2.*AL*AP*(tP( ? 1-1 Pill 1 ♦( AL2-AH2 MIM (3 I 
hCtDI ■hC)C-lB 1 


10 

zo 

30 

AO 

30 

GO 

.’0 

80 

90 

100 

no 

120 

130 

190 

150 

160 

170 

IBO 

193 

200 

210 

220 

230 

290 

250 

260 

270 

280 

290 

300 

310 

320 

330 

390 

350 

360 

370 

380 

39(1 

900 

910 

920 

930 

9 L 0 

950 

960 

970 

9S0 

990 

500 


IM 


CUM>OSm ANALYSIS ANL OESICN PRUiRAM > COHANC - C0HP4 


JULY# 19 7A 


2 5 T 

a 

o| 


ltl>NC1CT 

C LUNQlltUlKAL SINAIN CC'hSUAlSI - COtPatSSlON. 
CINCIOTi<LPl/lPLC>l. 
hCT01*KCID1*l 

C LONGlYLJlhAL STRAIN ClNSTRAINT - TENSION. 

CINCTCTI>EPl/LPLT-l. 

HCTCT«NCTLT#1 

C TRANSVLRST STRAIN CONSTRAINT - COnPRcSSIUN. 

6INCTCT>*lPZ/fcPTC-l. 

NCTCiI«NCTLT*l 

C TRANSyEkSL strain constraint - TENSION, 

CINCTLII«tP2/LPTT>l. 

KTOT>NCT(#Ta1 

C SHcAM strain constraint. 

GiNCTCTT'AbSUPJI/UNLT-l. 

IF (INNO.LT.N) GJ TO 160 
PAC'NAC 

CO 6C N*N1#NCTUT 
IF IGIP). Ut.cn NAC«PAC*1 
60 CUNTINLC 

IF (NAC.LG.NACI GO TO 160 
N2>NAC«1 
OJ 70 I1«N2#NAC 
CO 70 jj>;#NCv 
70 AI1I#JJI*C. 

NAC'NAC 

N2>N1 

DO ISC RN«1#NPIY 
r«Lnr(aki 

C GRACItNT CF <TmA|N> - LCUATION 37. 

00 60 N1>1#3 

TnpiNt i«u. 

OU «C A2«l#3 

BO THPCrtl I«1PPInH*CPli«l#N2#NA l♦c*’U2» 

CO 9C Nt«l#3 
OkPIRl I'L. 

00 90 a2-:#3 

JU Ol:P(Kn>OLPIKll-3dO(l#A2l*INPU2» 

CEPl>AL 2 * 0 tF{liAAR 2 *CEPI 2 lAAL«An« 0 cPI 31 
LcP2>AKjAlt P(ltAAk2*ulPI2l>AL*AnA0£P(3l 
0EP3#2.*A' *AM*IDEP(2l-0tPUnACAL2*AM2l«0EPI3l 
NaC’NAL 
N1«N2 

IF UlNll.LI.cn GO ID 100 

C CPAOIENT OF ACTIVE LDNCl T'JOINAL CONPRESSlVc STRAIN CONSTRAINT. 

NAL«NAC*1 

IF (NAC.EC.NN3I RETURN 
A<NAC#N i■All•AC•Al*OEPmPLC 
1C INAC •■Nl 
lOG Nl*Nl*l 


SIO 

520 

530 

5A0 

550 

560 

570 

5B0 

590 

600 

610 

620 

630 

6A0 

650 

660 

670 

6 B 0 

690 

700 

710 

720 

730 

7A0 

750 

760 

77C 

7B0 

790 

auo 

eio 

U20 
6 30 
DAO 
dSO 
660 
a 70 
860 
690 
900 
910 
920 
930 
9A0 
950 
960 
970 
980 
990 
luOO 


CJNPCSm ANALYSIS AND CtSlGN PROGRAtl - COnANO - COHPA 


JULY# 197A 



IF IGINll.lT.CT> CO TL 119 

1010 

C 

CRAOllNT Ct ACTIVE LUNulTUOlNAL TENSlLk STRAIN CONSTRAINT. 

1020 


NAC>NAC«1 

103- 


IF INAL.tC.NN31 return 

lOAO 


AINAC»N>«AlnAC#RI#UL*l/kPLT 

1050 


ICCNACI'N# 

1060 

110 

Nl»Ni*l 

1070 


IF (GINII. LT.cn Gb TO 12 J 

1080 

c 

GRADIENT b( ACTIVE TRANSVERSL COMPRESSIVE STRAIN CONSTRAINT. 

1090 


NAC>NaC*1 

1100 


IF INAC.cC.NN31 RETURN 

1110 


A(NAw#ll|«A«NAC#RI*uEr 2/EPTC 

1120 


ICINACl-Nl 

1130 

12 C 

Nl#Nl*t 

llAO 


IF ICINil.LI.cn GO TO 130 

1150 

c 

uRACUNT OF ACTIVE IRAUSVERSl UNSILE SIRAIN CONSTRAINT. 

1160 


NAC'NaC*! 

1170 


IF INAL.EQ.NH3t RFlUxN 

1180 


A(NAC#Nl •AlNAC>Rl*0EF2/tPTT 

1190 


IClNACl'Nl 

1200 

13C 

Nl.MAl 

1210 


Ir UtNll. lT.cn Gu T'. 150 

1220 

C 

gralieni cf aciive snfah strain CUNSTRAInI. 

1230 


NAC<NAE*1 

1250 


IF <NAC.Lb.NN31 RE*, URN 

1250 


S1GN>A. 

1260 


IF (EP3.LT.U.I SIGN#-’. 

1270 


A(NAC«l>l#AINAC>R|*SIbN>0EP3/GNLT 

12 B 0 


1C <NACI#n] 

1290 

lAO 

COMlNLc 

1300 

15C 

CONllNbt 

1310 

luO 

CONTINLE 

1320 

17C 

CDNTiNLr 

1330 

c 

cjnsipaints un stiffness. 

1350 


N1#NCIlT 

1350 


IF (AllL.LT.l.Jt-tuI GO 10 lorf 

1360 

c 

CONSIhAINT CN All. 

1370 


NC ICT#NC1CT«1 

1380 


CINCTr.TI#l.-AAIl.ll/AllL 

1390 

loC 

IF <A(2L.lT.1.9E-10> Cb TO 193 

1500 

C 

CJNSTFAIM LH A22. 

1510 


NCTLI#NCTCT#1 

1520 


G<NCTUTI>1.-AA12,21Va22L 

1530 

190 

IF (AclL.i1.1.9L-10I go IJ 2C0 

1550 

c 

LUNST#A|N| ON A66. 

1550 


NCTLI#NCTlT«1 

1560 


Gl KCTOT l#l .-AAI3# JI/At6L 

1570 

cOL 

IF <lNFU.lT.A.OR.M.tC.NCTJTI REIURN 

1580 


IF lAlU.ll.l.OE-lOl Cu T3 230 

1590 


Ni#M#l 

1500 


• 



, CLTHPOSITf AirALTSIS AND OESIDN PROGRAM - COnAHO - COHP4 JULr> 1974 



9 ^ 

S ^ 

S s 

.tl? 

g£3 

l| 

§ 


K1-NC1DT 

C LUNCIIDOJNAL STKAlh CDNETAAlNr - COHPRb:»StDN. 
CIKCTQT)>ePl/£PJ.C-l. ; . 

hCTDT-NCTbm 

C LOWDlTWINAt STRAIN CCN STB AtfiT - TCKSIOH. 

G(NCTDTl«Eni/LPLI-l. 

■ NCTOr-NCTUtl \ 

C; TRAKSVtRSI STRAIN CONS'] A AlnT - COnPRcSSIUN* 
ClNCJ0TJ-l.P2/t:PTi:-l. 

.. ! m:td]«nctdi+i 

. C : TRANSVEfiSe STRAIN CONSTRAlNr - tension. 

G(NCTtn«tP2/£PTTrl. 

NCTOT-NCTDTtl 

C SHEAR. STRAIN CONS TiIAIM;. 

C(NCTtT)>Afall£P3j>CNLT-iL* 

If tiNFO.LT.iVJ GO To 160 

. racAnac ■ • 

DO 6c K»M?NCTOT 

.■ ;lf (G{f:],6t,cn NACAPAC*! 

/ 6d CUNTliitE 

: IF (HACitQ.NAC} GO TO 160 

■ ■■ ha-NACti - ■ : 

DO 70 ll'NJfMAC 
' CO 70 jJ»i#NCV 
70 AtII,JJ]-C« 

MAC*NAC 
■: W2-M 

'DO 1st NNAlVhPLt 
: - K-LNRUKI ; 

C : GRADItNT CP STRATHS “ ECUATIOfl 37. 

DO 60 PlPlja 
TtiPt*ai-o, 

. DO ao KZ» 1 *D 

BO 7HP(KU.lPp<KH*CRlKl#K?>r.R)*ePtA2] 

DO 9C Rl«l,3 
DEPCKIJ.C. 

DO 90 RE-1^3 . : ■ 

90 0EPtltll*DEP(KU-33IKl#K2|*THPtrte> 

CEPl-Ata ♦ OEF tl HArti »CEP rai +ALFAfl*Q£? 1 3 ) 

. CtP2«AN2+LtPU}+AL(:*UfcPi2J-AU9AH*DEPC-3) 
DEPa*a,*ivnH*coEPtei“btpa])*-[ALa-ANai*DEPi3j 
NAC*HA(. ■ 

■ »il«h? 

: IF (GlNi). U1.cn Gii TO. TOP 

C CPADlEtlT Cp active LUNGtT'JOINAL COHPRESSlVc STRAIN CQNSTRAINr. 

:.NAC«NAC+1 : 

If tNAC.EC.NNSI RETURN 
, AihAC,fj*AU(ftc,K>+gep]/£Puc 
IC INACt-Nl 
lOG NI-NI+1 . 


510 
S20 
530 
640 
350 
660 : 
670. 
5B0 : 
590 
600 
610 
620 
630 
6^0 
650 
6S0 
670- 
6B0 
690 
700 
710 
7B0 
730 
740 
750 
760 
770 
7B0 
790 

eub 
BlQ 
U20.. 
B3D 
B40 
650 • 
660 
670. 
660 
690 
900 
910 
920 
930 
940 
950-: 
960 
970 
9BC 
990 
1000 
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IF (G1Ml.LT.CTi' CO TD 110 1010 

C GRADJLNtCf ACTIVE longitudinal TENSILL STRAIN constraint. . 1020 

NAC.-NACfl ICiaU 

. IF TNAU.tO.NH3): return-. . 1040 

AINAC»Kj-A{NAC.H>4UtPl/EPLT . ■ v 1050 

ICtNACi-M , TOGO 

.- 110 Nl-Nl+1 . ■ . .1071) 

IF ICIMJ.LT.CTl CD TO 120 1080 

C CRADlEMGf active TRANSVERStCDNORESStVE STRAIN constraint. - 1090 

NAC'NAC+I . : ■ Hoc 

IF INAC.tC.NN3l RETURN lllO . 

A(NAC»M'AUlAC»Rl*DET-2/tprC .. ■ 1120 

IC(NAC]>M .- 1130 

12C M«Nl+t 1140 . 

IP lOtND.LT.CII GO TO 130 1150 

C . oilADUNT Cf ACTlve rRAUSV.taSt TtHSILE STRAIN COHS.TRAENT. 1160 

NAC-Nacu 1170 . 

IF (MC.tO.Mr3t attURN. .. . IIBQ- 

A(NAC*M=iAtNACHt)>0EF£7tPTT 1190 

Iv(NAC1>n1 . 1200 ■ 

. ISO M^M+l 1210 - 

If fOir^D.tl.CTJ Gu Tf, 140 ; 1220 

C OHAtltNl CF ACTIVE .SNFAR STRAIN CJNSTRAIHT.' 1230 J’ 

NacANAC*!' 1240 ■ ■ 

IF fNAC.Ltt.MfS.] return . : ,1250 . 

storr'A. 1260 ' 

. IF .(EP3.UT.U. J SlGpi*-!. - 1270 

A(NACirPJ-A<NACfAJ+5ID[f*0C-P3/D6Lt . 12B0 

1C (NACJ.nl . 1290 

140 CUMlNUt -0; 1300 ' 

ISC CONUNUi : 1310 

lt/0 CilfiTlNLi 1320 

17C CUNTiNtf . ■ 1330 

C CPNSTPAINIS .liH SIIFFhESS, : : ; 1340 

Nl.NCIUT' 1350 

if [AilL.LTii.S'E''Iu) 00 tb'Ib6 . 1360 -. 

C CO NS THAI NT CN All, 1370 ; ^ ■ 

hCtOT.NCTCI*! " 1300 

GINtTDTI.l.-AAIl.II/AlIt 1390 :: 

loC. If 1A?2L.LT.1,0£-1D) GO TO 193 1400 : . ' 

C CUhSTfAtM CN A22. :,. 1410' 

NCUT.NCTCT+1 ' . 1420 

GtNCIU7).l.rAAi2f2J/A22L - 1430 

190. IF [Acc.L.ll.l.aE^101 Oj T-1 ZCO 1440 

C UONSTkAInI on A&6. 1450 

NCTuT«NCTuT+i ■ 1460 

U<RC;T0t>.I.-AS[3.sr/A66L - 1470 - ■ ; 

20C Jf (lijFU.LT.A.ilR.fil.IC.HCTJTj RETUaN 14 Bo 

If UllL.LT.l.OE-lOJ Gu TO 230 1490 

M.Nl+1 . . . 1500 ' .:■. 
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JULVf 197A 


IF IGlMi.LT.CTL) G'J 1U 210 

1910 

c 

SRAGIEKT C.F ACTIVl COHSIRAInT ON All. 

1520 

• 

naC>NACa1 

1530 


IF (NAl.tO.KKlt AtTURH 

19A0 


1C IKAC)*K1 

1550 


DO <10 K'lfKOV 

1560 

21C 

A(NAC>KI>C. 

1570 


DO 220 KR>i,KPiy . 

1580 


K'LNK (KK ) 

1590 

221) 

A(KAC>FI«A(KAC«l()-CR(lf 1*KRI/A11L 

loOO 

2J0 

IF (A22L.tT.1.0E-l0) GO 10 2oO 

1610 


Kl«hl«l 

1620 


IF (GIKll.LT.CTLt G1 TO 263 

1630 

C 

GRADIEM OF ACTIVE CCKSIRAtKT JK A22. 

1660 


KAC«KAC*1 

1650 


IF INAC.LG.KKII AhlURK 

1660 


ICINAC)>N1 

1670 


CO 2A0 «*l»KbV 

1680 

2<iG 

AtNACfRI«(. 

1690 


00 2tnv KK>lfKPLV 

'.700 


K'LKKIHKI 

1710 

2»0 

A(KAC»Pl>A(KAC»F)-CP(2f2»R<0/A22L 

1720 

26G 

IF lAttL.lT.I.OE-IO) RETURN 

1730 


Kl-M«l 

17A0 


IF (GtKlI.tT.CUl Gu TO 290 

1750 

C 

GRAOUKi OF AClIVt COK^IRAIKI JK A66. 

1760 


KAC>KAC*i 

1770 


IF (KAl.EU.KK3f RETlfRK 

1780 


IC (KACl'Nl 

1790 


DO 27C K*lfKDV 

1800 

270 

A(NAC>K)*(. 

IttlO 


DO 2eC hk'IiKPLY 

1820 


P<LKK(KR 1 

1630 

260 

Al KACf KI«A4KAC>H)-LFI3> If AKt /AObt 

1660 

zto 

COKTIKUE 

18 5C 


RElUftK 

I860 


END 

1870 
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Number 
of Cards 


COMAND DATA ORGANIZATION 




INFORMATION 


Title - Anything may be given here 


NCALC, NPLY, NDV, NLC, IPRINT 


LNK(I) ,I=1,NPLY 


X(I),I=1,NDV 


VLB(I) ,I=1,NDV (Blank card(s) if NCALC.NE.O) 


THN(I) ,I=1,NPLY 


EL, ET, GLT, PRLT 


EPLC, EPLT, EPTC, EPTT, GMLT, AllL, A22L, A66L 


PN(J, I) ,J=1, 3 (One card per loading condition) 


Begin with next set of data - Program terminates 
if 2 blank cards are read here. 


FORMAT 



8F10.2 


8F10.2 


8F10.2 


4F10.2 


8F10.2 


3F10.2 



TABLE 1 - DATA ORGANIZATION 




































BiCTA SHEET 


COMPOSITE ANALYSIS AND DESIGN PROGRAM COMAND 



1 5 

11 

21 

“1 

31 

41 

IHi 

61 

71 

A 






B 


— 










C 



■■ 













• 











D 



























E 



























P 



















■ 








G 







■ ■ 


H 









I 


































. 

1 












TABLE 1 - DATA ORGANIZATION - CONCLUDED 
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JANUARY 1973 


TABLE 1 . 2 . 1 -in. KEY UNIDIRECTIONAL PROPERTIES 

I HtGH-STRENGTH GRAPHITE/EPOXY - [0) 


Vf = 0. 60 



RT 

350*F 

Design 

strengths* 

Longitudinal tensile ultimate 

Ksi 


180. 0 

180. 0 


Transverse te:»sile ultimate 

Ksi 


8. 0 

4. 0 


Longitudinal compression 
ultimate 

Ksi 

^1" 

180. 0 

70. 0 


Transverse compression 
ultimate 

Ksi 

T 

30. 0 

12. 0 


In-plane shear ultimate 

Ksi 

F^u 

LT 

12. 0 

6.8 


Interlaminar shear ultimate 

Ksi 

j,isu 

n. 0 

8. 0 


Ultimate longitudinal strain 

flin. /in. 

tu 

‘l 

8, 700. 0 

9,650. 0 


Ultimate transverse strain 

fiin, /in. 

‘■t" 

4,750.0 

4, 100. 0 

Elastic 

properties 

[typical] 

Longitudinal tension modulus 

Msi 

< 

21. 0 

18.7 

Transverse tension modulus 

Msi 

4 

1.7 

0.87 


Longitudinal compression 
modulus 

Msi 

^L 

21. 0 

18. 7 


Transverse compression 
modulus 

Msi 


1.7 

0. 87 


In-plane shear modulus 

Msi 

‘^LT 

0. 65 

0. 32 


Longitudinal Poisson' s ratio 



0. 21 

0. 21 


Transverse Poisson' s ratio 



0. 017 

0. 010 

Physical 

Density 

lb /in. ^ 

P 

0. 056 

0. 056 

constants 

[typical] 

Longitudinal coefficient of 
thermal expansion 

^in. /in. /"F 

«L 

-0, 21 

-0. 005 


Transverse coefficient of 
thermal expansion 

//in. /in. /*F 

“t 

16. 0 

21. 8 




.-i 




References: 1.2-15, -19, -21 

^Typical Design Allowable, reference section 1.2.0 

1.2.1 Table 2.- Material properties. 
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G 
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.21 





H 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

I 

180000. 

0. 

0. 







0. 

-30000. 

0. 








8000. 

0. 







0. 

0. 

12000. 
















TABLE 3 - DETERMINATION OF LIMIT STRAINS - G/E COMPOSITE 
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TABLE 5 - (0, ±45, 90) GRAPHITE EPOXY COMPOSITE - EXAMPLE 2 

















































DATA SHEET 


COMPOSITE ANALYSIS AND DESIGN PROGRAM - COMAND 


'“'---..Col 
Block'“ — ^ 

1 5 

11 

21 

31 

41 

51 

61 

71 

A 

CO. 30, -3C 

1. 60. -60. 

90) GRAPH] 

[TE EPOXY r 

OMPOSTTE - 

EXAMPLE 3 . - 

B 

0 

6 

4 

4 

0 







C 

1 

2 

2 

3 

3 

4 



1 

1 - 







— 
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.1 

.1 

.1 

.1 
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.00001 

.00001 

.00001 

.00001 
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0. 

30. 


60. 

-60. 

90. 
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,21 
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.00471 
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TABLE 6 - (6, ±30, ±60, 90) GRAPHITE EPOXY COMPOSITE - EXAMPLE 3 
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Figure 1.- Typical composite loading. 
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Figure 2.- Typical ply orientation. 
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Figure 3.“ Determination of limit strains - G/E composite 
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Figure 4.” Quasi- isotropic G/E composite under 
uniaxial load — Example 1. 
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Figure 6.- (0, ±30, ±60, 90) Graphite epoxy composite — Example 3. 
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gu-e 8.- Ultimate compressive strength high-strength graphite/epoxy - 


°i/±**5j/9oJ family. 





su 

Allowable strength F - Ksi 


47 


(- 


80 

70 

60 

50 

ko 

30 

20 

10 


0 20 40 60 80 100 120 

i ±45® Plies 

Figure 9*“ Ultimate shear strength F^ high-strength graphite /epoxy - 

Ajr 

(0i/±U5,)/90k| family. 
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FIGURE 11.- 'COMAND' BLOCK DIAGRAM. 








